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Abstract 

A derivation of the Cesaro-Fedorov relation from the Selberg trace formula 
on an orbifolded 2-sphere is elaborated and extended to higher dimensions using 
the known heat-kernel coefficients for manifolds with piecewise-linear boundaries. 
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1. Introduction. 

In an earlier work [1] we have indicated how the relation between the numbers 
of symmetry planes and symmetry axes of a solid discovered by Cesaro [2], and 
Fedorov [3,4] can be derived analytically. In the present work we wish to further 
analyse the method and extend it to higher dimensions. 

This work can be looked upon as a discussion of the heat-kernel and (function 
on certain manifolds and so should be relevant to quantum theory. 

Cesaro's proof of the formula relies on an earlier result of his [5] relating the 
number of symmetry planes to the number and order of symmetry axes lying in a 
single symmetry plane. For completeness we repeat his findings here. Let b\ be 
the total number of symmetry planes and let there be, on the i-th such plane, rim 
symmetry axes of order to, then Cesaro [5] shows that (see later) 

&i = l + 5>g?(m-l), Vt. (1) 

Summing this formula over all symmetry planes and using the fact that 

J^rajJ = mn m , 

i 

where n m is the total number of symmetry axes of order to, we find one form of the 
C-F relation, 

b 1 (b 1 -l) = J2 n mm(m-l). (2) 

m 

Fedorov [3,4] refers, somewhat confusingly, to x(x — 1) as the 'function' of x. (2) is 
his 'zonohedral formula'. (See [6].) 

Another form of the relation follows by noting that, for the regular solids at 
least, we have the orbit-stabiliser result 

g = 2mn m , Vto, (3) 

where 2g is the order of the complete symmetry group, V. Therefore 

2h (&! - 1) = g ^(m -l)=g(p + q + r-3) (4) 

m 

for the three types of axis of orders p, q and r. 

Relations (2) and (4), and their generalisations, are our interest in the present 
paper. Because the theory of polytopes might not be too familiar, a certain amount 
of historical and technical commentary might be of benefit, if only to the writer. 
Everything goes back to Schlafli. 
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2. A survey of the geometrical method 

We recall the standard geometric situation [7], [6] §3-5, §4-5. The regular polyhe- 
dron {p, q} is projected onto its circumsphere giving a spherical tessellation. The 
planes of symmetry of {p, q} divide this tessellation further into |T| congruent spher- 
ical (Mobius) triangles, (pqr) (with interior angles ir/p, ir/q and n/r), which are 
transitively permuted under the action of T, = [p, q]. This produces a simplicial 
decomposition (triangulation) of the sphere. The 'g- vertex' (0) of (pqr), and its 
images under T, correspond to the vertices of {p, q} and the 'p-vertex' (2) likewise 
to the vertices of the reciprocal solid, {q, p}, i.e. to the mid-points of the faces of 
{p, q}. The 'r-vertex' (1) is associated with the mid-points of the edges and, for 
the usual solids, r = 2 so that g = 4n r . The edge mid-points form the vertices of 

the quasi-regular polyhedron j P j which projects to a spherical tessellation of edge 

length 2ir/h, where h is the Coxeter number ([6] §2-3, §2-6, 5-91). 

Equation (3) says that there are g/m images of the m-fold vertex of (pqr) under 
[p, q] . This follows because the stability group of this vertex has order 2m. 

Steinberg [8] gives a simple geometrical derivation of (2) which appears in the 
second edition of Coxeter's book [9]. In the first edition, [6], Coxeter just indicates 
the general nature of such a proof. Curiously, in [9], he drops the references to 
Cesaro and Fedorov, presumably because their relations are superseded by Stein- 
berg's and some further work of his own, [10]. 

The geometrical argument is worth repeating. The b\ symmetry circles inter- 
sect in b\(b\ — 1) (possibly coincident) points (the vertices of the triangulation). 
In order that a vertex correspond to a symmetry axis of order m, it is necessary 
that m symmetry circles cross there. That is to say, there must be a coincidence 
of m(m — l)/2 intersections at the vertex. Equation (2) expresses the distribution 
of the total number of intersections among such vertices, and there are two vertices 
per axis. We might term m(m — l)/2 the degeneracy of the intersection. 

It is clear that this argument is essentially the same as Cesaro's. Indeed (1) 
follows from a counting argument applied to a single reflecting circle. The total 
number of such circles is just one more than the number that cross a chosen circle. 
Further, trivially, the number that cross this circle at a vertex of order m is m — 1 
(discounting the original one). 

We might note the relations between the n m and the numbers of vertices, edges, 
and faces, N , N± and N 2 respectively, of the polyhedron {p, q}, ([6] 5-81), 

2n p = N , 2n q = N 2 , 2n r = N ± . (5) 

The statement that qN = 2Ni = pN 2 is a graph-theoretic one (see e.g. [11] 
Chap.8). 

It is worth outlining the further development for completeness [6], [8]. The 
associated polyhedron |^| has Ni vertices, 2Ni edges (from graph theory) and 
Nq + N 2 faces (of two sorts, {q} and {p}). Every edge belongs to just one equatorial 



3 



/i-gon of which, therefore, there are 2N±/h in total. Since there are four faces at 
each vertex, the intersections of the /i-gons are simple (degeneracy of one) and the 
total number of intersections is just the total number of vertices, i.e. 

Equivalently, the number of equators is one plus the number of intersections with 
one equator i.e. 2N\/h = 1 + h/2. Thus 

2g = h(h + 2). (6) 

It is easily seen that each equator on {p, q} contains h points 1 and crosses h 
segments 02. The total number of reflecting circles on {p, q} may now be obtained 
by adding up the crossings on a complete circumnavigation of an equator. There 
are two for each 1 point and one for each 02 segment crossed. Taking into account 
the antipodal points we get b\ = (2h + h)/2, 

From (7), (6) and (4) we find Steinberg's formula [8], [9], 

24 

h + 2 =T7 ) (8) 

12 — p — q — r 

for the Coxeter number in terms of the angles of the Mobius triangle. 
3. The analytical method 

Let ((s) be the (^-function for the Laplacian on the orbifolded 2-sphere, i.e. on the 
Mobius triangle (pqr). The analytical method consists of identifying £(0) calculated 
in two ways - one from the expression in terms of the degrees d of the symmetry 
group, r, (this involves the Barnes (^-function) and another that comes either by 
using the expression for the constant term in the heat-kernel, derived by differential 
geometric methods for any manifold and then applied to the orbifolded sphere, or 
by using the particular form of the (^-function on the orbifolded sphere in terms 
of the symmetry axes, i.e. fixed points, as given in [12]. We could say that the 
Cesaro-Fedorov (C-F) formula is an application of the Selberg trace formula to these 
spherical domains. 

The formula referred to in [12] gives a constructional form to the (^-function in 
terms of 'cyclic' (^-functions. As derived, the formula is for the rotational subgroup 
r + and a special conformal coupling. It comes originally from the conjugacy class 
decomposition of the heat-kernel and is actually independent of the coupling. 

Although convenient, using this particular form of £ on S 2 /r is, in some ways, 
not as satisfying as an appeal to a more universal formula such as that for the 
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constant term in the heat-kernel expansion used in [1]. We need no knowledge of 
group theory nor of conjugacy class decomposition to use or derive this formula. Of 
course, the significance of the corner terms (see later) is, at root, a fixed point one. 

The number of reflecting planes is introduced analytically through the results 
of Shephard and Todd [13] and Solomon [14] who give this number as the sum of 
the exponents, (di — 1), of the finite reflection group. Our approach should thus be 
valid for solids other than the regular ones. 

Clearly the C-F relation would emerge whatever the coupling (we shall verify 
this shortly) and also there is no need to use the full, reflecting V. We could have 
obtained the C-F formula from the results of [12]. 

Another point we shall investigate concerns the value of the (^-function at other 
arguments, in particular at negative integers. One would expect to get different 
relations between the orders m and the number of symmetry planes. 

4. Conformal— rotational derivation 

In this section we rederive the C-F relation using the special conformal coupling 
£ = 1/8 and the corresponding rotational ^-function on S 2 /r derived in [12]. We 
refer to equation (26) of [12] which we give again with a slightly different notation 



CM - \ 



^]mn ra ( m (s) -2( J^n m - 1)Cjj(2s - 1,-) 



(9) 



Cm(s) is the (^-function for the rotational cyclic group, Z m , and its values at 
the negative integers are shown in equation (32) of [12]. (See the appendix to the 
present paper). In particular we calculate Cm(0) 



Cm(0) = 



2m 2 - 1 
12m 



(10) 



and so we get from (9), 



C(o) = 



12<7 



2^n m (m 2 -1) + 1 



(11) 



On the other hand we know that the rotational £(s) is given by the sum of the 
Neumann and Dirichlet (^-functions, 



where 



((s) = (n(s) + (d(s) 
d-l 



( N {s) = ( d (2s, 



d) 



( D (s) = Q(2s,Zd t -^-±\d). 



(12) 

(13) 
(14) 
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For s a negative integer, or zero, the two terms in (12) are equal. 

Barnes [15] gives the values, (we use both Barnes' notation and the more 
standard one of Erdelyi, [16], for the generalised Bernoulli polynomials) 

«-*■» i d > = t+t = Tpr(^ <)t( ° 1 d) (16) 



where 



w a~ L v-/ -1 V^- ~7 211* ' 

(a) = ^; d - 1} (£d, - a) = (6a 2 - 6a £ d, + £ d 2 + 3 £ d.dj) 

(16) 

and so we get from (12), (15) and (16), with g = d\di and b\ = d\ + d 2 — 1, 

C(0) = ^ [3 - 6(dx + d 2 ) + 2(d 2 + d 2 + 3^)] = ^- [2 0l ( 0l - 1) - 1 + 2o] . 

(17) 

The equation that must hold, therefore, is 

b 1 (b 1 - 1) = 5> m (m 2 - 1) + 1 - g. (18) 

Recalling the total rotational order equation 

Y J {m-l)n m = g-l (19) 

m 

we get 

&i(6i - 1) = ^n m m(m - 1) 

m 

as required. 

Compared with Steinberg's geometrical argument, this is a rather complicated 
derivation. Actually, it is not so very different in content. We are simply re- 
extracting the information that has already gone into the construction of the £- 
function (9). This remark leads onto the next point. We note that the C-F relation 
has been derived in its first form, (2). The derivation in [1] yields the second form 
(4). Let us see why this is so. 
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5. Heat kernel method 

For the differential operator —A 2 + ^R, the general formula for £(0) on a two- 
dimensional domain, .M, with boundary dM. = UdMi is 

«°> ^L r "^w^ L K(i) * + ^ ? ^ (20) 

where the # sum runs over all inward facing angles at the corners of dM. For 
minimal coupling £ = 0, while for the special conformal coupling in [12], £ = 1/8. 

The corner term was derived by Brownell [17] and Fedosov [18] and again by 
Kac [19] as the constant term in the short-time expansion of the heat-kernel. The 
angles 9 are quite general. 

Restricting attention to the spherical triangle, M = (pqr) , k must vanish since 
the sides are geodesic. Further, if the triangle tiles the two-sphere once, 

^ t RdA — 1 - ^ 



24tt J m (>g 

2g being the number of triangles. 

Choosing £ = 1/8 (not essential), so that £(s) = Cn(s), the identity becomes 



e 

or, writing the angle 9 typically as jv/m 



(21) 



*£( m -^) =2(&i(& 1 -l)+s-l). 



(22) 



Reorganising the left-hand side slightly just to see how things work out in detail 
g J2(m -l) + Sg- g Y t - = 2(b 1 (b 1 -l) + g-l). (23) 



m 

m 



The Gauss-Bonnet, or area, equation reads 



m q 



and we see that the 2(g — 1) term on the right of (23) will cancel against a term on 
the left to regain (4), bypassing explicit mention of the numbers n m . We do have 
to know the area formula for a spherical triangle, however. 
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The n m can be reinstated using (3), when we see that (19) is equivalent to 

(24). 

For future reference, we note that the quantity ^ — 1 is the number of elements 
of the complete symmetry group that fix an axis but not a plane in R 3 . It is just 
the number of nontrivial rotational elements i.e. all the elements in F + except the 
identity. We denote it by b 2 and write (19) as 

b 2 = J2 n m(m-l) (25) 



which can be looked upon as an identity similar to, but more elementary than, (2). 

Further numerical expressions, which are easily obtained, give the numbers of 
vertices, V, edges, E, and faces, F, of the triangulation in terms of the group order, 
(cf [6] p 67) 

V = 2 + g, E = 3g, F = 2g. (26) 

Equation (24) can be interpreted in the following way. Think of the combina- 
tion 2g 5^(1 /2m) as a sum over all 2g fundamental domains of a quantity, ^(l/2m), 
associated with one domain. If 2g Mobius triangles are assembled to form the sim- 
plicial decomposition of S 2 , each vertex will appear with a multiplicity of 2m. The 
l/2m factor reduces this multiplicity to unity and so 2<7^(l/2m) equals the number 
of vertices, V. 

6. Other relations. Higher dimensions 

Other relations can be derived using (12) and our general knowledge of the heat- 
kernel coefficients, or the specific form of the (^-function on S d /r, such as (9). 
The standard facts we need are that, if the heat-kernel expansion is 

K W*7&m E ^ < 27 > 

v ' fc=0,l/2,... 

then 

= ( 4n y/2 Ck + d / 2 ~ n o$ko, (28) 

where no is the number of zero modes, and, further, that Cd( s ) has poles at s = 
(d — m)/2, for m = 0, 1, . . . , d — 1 and m = d + 1, d + 3, . . ., with residues 

Un)- d / 2 

-C m/2 . (29) 



r((d-m)/2) 

These equations can be read in the two obvious ways - knowing the values of 
the ^function we can derive the coefficients and vice versa. For the moment we 
adopt the latter viewpoint and use the expressions for the coefficients developed 
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using standard analytical methods. Our discussion is neither comprehensive nor 
systematic. 

Two basic dimension- and coupling-independent results are 



Co 



\m\, c 1/2 = ±y^j2\ dM ^- ( 3 °) 



The ± signs are for Neumann and Dirichlet conditions. 

From (29) we see that Co and Ci/ 2 are given by the residues in £(s) at s = d/2 
and s = (d — l)/2 respectively. For simplicity we might as well use the conformal 
C-functions given in (12). Therefore we need the pole structure of the Barnes 
function [15] which is, as s — > k/2, 

C " (2 "- a|d) ^ 2(fc-l)! s-k/2 (31) 

where k is one of 1, 2, 3, ... , d. 

Comparing (31) and (29), cf [12] eqn. (76), 

C m/2 = s [ d - m +1) (a) 

7 V ; r((d-m+l)/2) 1 V ; 

where a is (d— l)/2 for Neumann conditions and ^ di — (d — l)/2 for Dirichlet. For 
simplicity, we usually work with Neumann conditions. In view of (16) and (30) we 
get in two dimensions, for example, 

\M\ = ^r, \dM\ = 2n dl+ /r 1 = 2 -^- 03) 
d±d 2 did 2 g 

These equalities are easily confirmed from the trivial geometrical fact that the 
two-sphere is covered by 2g spherical triangles, the vertices of which are formed by 
the intersections of the b± symmetry great circles. The sum of the circumferences 
of these circles is half the total perimeter of 2g fundamental domains, as is quickly 
appreciated by drawing in the sides of each triangle separately - equivalent to 
splitting each great circle into two (cf [6] p232). More formally, one can use the 
measure-theoretic principle 



\AUB\ + \AnB\ = \A\ + \B\ (34) 



applied to the edges. 
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Dividing in (33) we find 

\BM\ 

h = W 

which is Laporte's rule, [20], for the number of symmetry circles, i.e. the nodal 
lines of the lowest odd mode (the Jacobian of degree b\). 

The extension of these elementary results to higher dimensions is immediate. 
The relevant analytical fact is that, for all d, 

d S[ d+1) (a) = -, d s[ d) (E^-(d-l)/2) = ^ (35) 
9 ^9 

so that 

\M\ = -^-\S d \ (36) 

and 

\dM\ = —\S d - 1 \. (37) 
9 

The generalised Laporte rule is then 

|S d | \dM\ 

&1= Jcd-u 7771 • (38) 



The geometric situation in is the generalisation of that in M 3 , recalled in 
section 1. We outline it here for any d (see [21], [6] pp 130, 137-140, [22] §4-1, [23] 
chap. 12). We are most interested in d = 3. 

The projection of the (d+ 1) -dimensional regular polytope, {ki, &2, . . . , kd} (in 
IR d+1 ) onto its circumscribing hypersphere yields a spherical tessellation, or hon- 
eycomb, the cells of which are the projections of the d-faces of the polytope. The 
symmetry d-ft&ts ('reflecting hyperplanes') of {k\, &2, . . . , kd} intersect the circum- 
sphere, S d , in a set of reflecting great (d— l)-spheres constituting the boundaries of 
|T| spherical (i-simplexes that are transitively permuted by the complete symmetry 
group T, = [ki, /c2, • • • , kd] and which make up a simplicial subdivision of S d , [24]. 
(This decomposition could also be termed an irregular tessellation.) The (d — 1)- 
dimensional boundary of the fundamental domain is the union of d pieces of the 
reflecting great (d — l)-spheres. Topologically, the honeycomb is the boundary, or 
frontier, of the 'solid' polytope. For d = 3 the fundamental domain is a spherical 
tetrahedron. Goursat [25] gives a detailed account of the different S 3 honeycombs. 

Note that the symbol {k±, fo, • • • kd} can refer to a (<i+l)-dimensional polytope 
in M d+1 or to a spherical (d + l)-dimensional polytope in S d+1 or to a (i-dimensional 
honeycomb on S d , depending on context ([6] p 138). 

Equation (37) expresses the fact that the measure of all the reflecting great 
spheres is one half the measure of the total boundary of 2g fundamental domains. 
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The next identity arises on setting m = 2 in (32), 

d = 2n\S d - 2 \ d S i 1 d ~ 1) ((d-l)/2). (39) 

The C\ coefficient in two dimensions has already been discussed in section 2 
but we wish to extend the argument to <i-dimensions and so write out its general 
form, valid for any coupling and region M., 

Ci= 1 —^ I RdV+\j2f ndS+\y^' K -—^\dM l ndM J \. (40) 

b JM 6 i JdMi b ~T* Vij 

Here 6^ is the (arbitrary) dihedral angle between the boundary components dM.i 
and dAij. For simplicity we have assumed that this angle is constant along the 
intersection. This is so for a fundamental domain. Fedosov's calculation, [18], of 
the final term is for a flat polyhedral region so that the 9^ are constant for him. 
For d = 3 the term could be referred to as an edge term. 

We again employ the conformal coupling of [12] so that £ = (d — l)/4d. For 
M. a fundamental domain, the first term simplifies to 

1^1 jMK-C-'W-Vl (41) 

and, because the boundary is geodesic, the second term in (40) vanishes. 

We now look at the right-hand side of (39) which can be written more geomet- 
rically in terms of the number of reflecting (i-planes and the number of those group 
elements that fix a (d — l)-plane but not a <i-plane. We denote these numbers by 
bi and 62 respectively. A reflecting (d — l)-plane intersects the circumsphere in a 
reflecting great (d — 2)-sphere. 

The complete symmetry group V is generated by (d+1) reflections, Ri, . . . , 
Rd+i, in (d+ 1) hyperplanes in Generally, b r is the number of elements of F 

that fix a (d + 1 — r)-flat, but no flat of higher dimension, in IR d+1 . Equivalently, 
b r is the number of elements that are expressible as products of r (but no fewer) 
reflections, i.e. those that have length r. 

The Coxeter element, R1R2 ■ ■ .Rd+i, has period h, the Coxeter number, and 
characteristic roots exp(27rimj//i) where the rrij (j = 1, . . . , d+1) are the exponents, 
related to the degrees by rrij = dj — 1. 

Solomon's theorem, [14], states that 

d+1 d+1 

\[(l + m i t) = Y J Kt r . (42) 

i=l r=0 

We have chosen the last exponent, m^+i, to be the one that equals unity. Therefore 
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d+1 d 



d d 

i<j 



b 2 = ^ miirij + ^ ' 



and 



d+l d+l 
63 = mimjmk, 64 = mimjirikmi. 

i<j<k i<j<k<l 

Because the characteristic roots occur in complex conjugate pairs, we have 

61 = J2 mi = (d+l)h/2, (43) 

in terms of the Coxeter number (cf (7)). This result also follows from a generali- 
sation of the equator argument ([9] p 229-231, [26]). We also note that the largest 
exponent equals h — 1. 
The breakdown 

2g = J2br 



r 



is obvious. Setting t = — 1 in (42) we have 

£(-l) r & P = (44) 


and so 

9=J2 b r = J2 b r- ( 45 ) 
even odd 

The first sum gives the order of the rotation subgroup, and the second the (equal) 
number of remaining elements. 

Returning to (39), simple algebra yields 

((d - l)/2) = ±- (hib, - 1) + b 2 + (3 - d)/2). (46) 

(The systematic evaluation of the Bernoulli polynomials is discussed in the ap- 
pendix.) 

When (46) is substituted into (39), the last term cancels against (41) finally 
giving the identity 



if A V ij 77 / 



(47) 
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As a simple check we can set d = 2, when b 2 = g — 1. Using |S°| = 2, (47) 
reduces to (22). (As a minor point, we note that for d = 2, the identity does not 
come from residues.) 

For d = 3, the easiest case to visualise, the dAii, (i = 0, 1, 2, 3,), are the 
spherical triangular faces of Ai, the spherical tetrahedron fundamental domain of 
the simplicial decomposition of the 3-honeycomb II4, = {k±, k 2 , k^}. The edges, 
dAii fl dAij, are the six, circular sides of these faces. The situation is depicted in 
[6] Fig. 7-9A, [22] Fig. 4-3B, [25]. Verbal descriptions can be found in [9] p 229, 
Todd [27] p 216, and Sommerville [21] p 188, [28]. 

The vertices Pi (i = 0, 1, 2, 3) of the fundamental domain are obtained as 
follows. P3 is the centre of a cell of the honeycomb, i.e. of the projection onto the 
circumsphere of a 3-face of the 4-polytope, II4. This cell is itself a regular (spherical) 
polytope, II3, and P2 is the centre of a cell of its boundary. This cell is a 2-face 
('face') of II3 (and of II4) and is also a regular polytope, II2. The cells, Ill's, of its 
boundary are 1-faces ('edges') and Pi is the centre of one of them. Finally, P is 
the end of such an edge and is the one point that corresponds to a vertex of the 
original polytope, 1I4. An edge of a cell of the honeycomb has length 2|PoPi|. A 
nested set of TI^ (i = 0, 1, 2, 3) is called a flag. Cells are sometimes referred to as 
facets. 

Goursat [25] gives perhaps the most comprehensive description of the geometry. 
He refers to elements of the rotational subgroup, r + , as 'transformations droites' 
(i.e. even) and to the remaining elements as 'transformations gaudies' (i.e. odd) 
and expresses them in terms of, up to four, inversions. 

If a face of Ai be labelled dually by the opposing vertex, then the dihedral 
angles, between faces i and j are given by #01 = ir/ki, Q\ 2 = n/k2, 6*23 — n/ks 
with the rest being 7r/2's ([27], [6], [28]). 6ij can also be characterised as the 
dihedral angle opposite to the edge PiPj- 

We call the edge lengths L^, = |PjPj|, and expand the identity (47) for d = 3 

as 



To check (47), or (48), we need the edge lengths. Some of these are given in 
[6], p 139 and Table I (ii). 

It is always the case that Lqi + £23 + -^03 an d £12 + -^02 + -^13 are separately 
commensurable with n. For example, in {3 3 } these quantities are both equal to it 
(so that the total perimeter is 2tt) and (48) is then easily checked numerically using 
g = 60, bi = 10 and 62 = 35. 

General formulae for L^- exist due to Sommerville [28]. They are handy nu- 
merically and, for convenience, are repeated here in the appendix. We have used 
them to confirm (48) in all cases. 




(48) 
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7. Geometrical meaning 

A blind, numerical verification is all very well but equation (47) has a geometrical 
interpretation, or, if one prefers, a geometrical derivation. The left-hand side is 
clearly a sum over all the fundamental domains of a quantity associated with one 
domain and the right-hand side must be a global measure of the same total quantity. 
We know, for example, that 61(61 — l)/2 is the number of S d ~ 2 intersections of b\ 
S d_1, s, including coincidences. The corresponding identity will be derived from (47) 
to see how the various terms fit together. Our treatment is not the most direct. 

For simplicity, attention is restricted to d = 3. It is convenient to depart from 
the traditional notation and define 

7T 

k k i = —, (ijkl) cyclic. (49) 
Vij 

The two-dimensional case suggests that we seek to evaluate the term 

-Lki 



in (47). Assembling 2g fundamental domains, the l/2kki factor cancels the multi- 
plicity of appearance of the arc PfcPz and the result is the total circumference of all 
the circular intersections of the reflecting S 2, s. Hence 

2nb' 2 = gJ2l^ L ki (50) 

k<i Kkl 

where b' 2 is the total number of complete S 1 boundaries in the simplicial decompo- 
sition of the honeycomb. [For {3 3 }, b' 2 =25.] 

The same reasoning applied to the perimeter gives the formula 

2gY J L kl = 2Y J N kl k kl L kl (51) 

k<l k<l 

where N kl is the number of PfcP; arcs in the simplicial decomposition. This agrees 
with 

N kl = ^- (52) 
kkl 

which follows directly from the orbit-stabiliser theorem. Using (50), (47) becomes 
(for d = 3) 

g k kiLu = 2tt(6i(6i - 1) + b 2 + b' 2 ) . (53) 
k<i 
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Define, for convenience, 

D = gJ2Lki-2n(b 2 + b' 2 ) (54) 

k<i 

to get the intermediate result 

2nb 1 (b 1 -l)=gJ2( k ki-l)Lki + D. (55) 

k<i 

We will show that D — 0, proving, incidentally, that the perimeter of the fundamen- 
tal domain is commensurate with n. We already have this because the perimeter is 
half the sum of the four face-perimeters each of which is known to be commensurate 
with 7r. We then see that b 2 + b' 2 is the number of S^s on all the bounding S 2 's, 
counting as distinct, circles that belong to two, or more, S 2, s. 

Proceeding in the geometrical vein, and noting that the intersection degeneracy 
of the reflecting S 2, s at the (kl) -edge is kki {kkl — l)/2, the total perimeter of the S 1 
intersections, allowing for coincidences, is 

2tt& 1 (& 1 - 1) = Y, Nklk ki{hi ~ l)L kl . (56) 

k<i 

Hence, using (52), a comparison with (55) shows that D = 0, which equality can 
be put into the form 

2nb 2 = J2 NH (kki-l)L kl . (57) 

k<i 

We have thus obtained identities, (56) and (57), precisely analogous to (2) and 

(25). 

In M 4 the analogue of a rotation axis in IR 3 is a rotation plane, i.e. an S 1 
boundary on the circumsphere, S 3 . Such a circle is divided into a number of arcs 
each with an associated order of rotation, one of the k k i- In words, (57) says that 
62 is obtained by counting the number of nontrivial rotations at every arc, weighted 
with that arc's relative size. 

Although we can look upon our development as a proof of (57) it ought to be 
derivable directly. However, we have been unable, so far, to provide a geometrical 
reason for (57). It might be related to the various extensions of the Euler formula 
such as the Dehn-Sommerville equations and other angle-sum relations, of which 
there seem to be quite a variety, e.g. [29]. 

To set the geometrical picture more closely, the general situation is again re- 
called. 

Each reflecting (d— l)-sphere is divided by its intersections with all the others 
into a certain number of (d— l)-simplexes taken from those forming the boundary of 
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the fundamental domain. The topological structure is the same for all the reflecting 
spheres although the composition may differ metrically. 

From (37) we find that the number of (d — l)-simplexes in one reflecting (d — 1)- 
sphere is (d + l)g/b\. If we assume (43), this number equals 2g/h, ([6], [9] pp 
231-232, [26]). 

A simple example is the {3,4} 2-honeycomb for which there are two types of 
reflecting circles. One has the vertex content 02120212 and the other 01010101. 
(See [6] p 65.) Both have 2g/h = h + 2 = 8 segments. Incidentally this is a good 
example in which to check (1). 

We need to prescribe the polytope in more detail. This is done in terms of 
its configurational numbers A% ([6] pp 12, 130, [21], [28], [23] §12 • 6). A regular 
polytope, Il n can be thought of as a nested sequence of regular polytopes, of type 
11;, (i = 0, 1,..., n). (See [6] pp 127-129.) An element, Uj, belongs to N jk II fe 's 
for each k > j and contains Njk life's for each k < j. The number iVj, which for 
symmetry's sake can be written N ii: is just the number of Ilj's in the polytope. N 
is the number of vertices, for example, and N n is one. 

For example, the regular simplex, aj+i = {3 d }, has a complete symmetry 
group of order (d+ 2)! and numbers, N, 



J + 

(See [21] p 96, [6].) Notationally, Coxeter in [30] has ( l \d+ l) for N t . 

If we set n = d + 1 and take Hd+i to be the Euclidean polytope that projects to 
the (i-dimensional honeycomb, then the vertices, Oi, (i = 0, 1, . . . , d) of the (d+ 1)- 
simplex in the decomposition of Hd+i project to the vertices Pi mentioned earlier 
(for d = 3). Od+i is the centre of Hd+i- 

For the (d + l)-polytope Hd+i = {ki, ki-, ■ ■ ■ kd} we can put 

Ui = {fci, k 2 , • • • , 

or, for the ascending construction, 

II- = {kd-i+2, kd-i+3, • • • , k d }, 

for the constituent polytope elements. The relation between the orders of the com- 
plete symmetry groups of these IIVs and the numbers N is 

1 \[h, . . . , ki^ll . \[k i+2 , . . . , kd}\ [oy> 

which follows from the orbit-stabiliser theorem. The first factor in the denominator 
is the order of the 'internal' symmetry group of 11^ while the second is the order of 
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the stabiliser of Ilj as a whole. The ratio is then the index of the subgroup leaving 
the centre of Hi, that is P^, or Oj, invariant i.e. it is the number of n^'s. We must 
remember that the internal symmetry group of an edge, {}, has order 2. Equation 
(58) is an example of (59). We note that Ni = N' d _ i . 

Two special elements are 11^, which is the 'cell' of lid+i (sometimes called the 
'bounding figure'), and 11^, which is the 'vertex figure'. A property of the regular 
polytope is that the vertex figure of the bounding figure is the bounding figure 
of the vertex figure. Continuing with the terminology, we may refer to 11^ as the 
(d + 1 — i)-th vertex figure (cf [6] pp 133, 134.) and the Ilj as the (d + 1 — i)-th 
bounding figure. The 11^ are the perischemons of Schlafli. 

8. Volume and degrees 

A fundamental result of Schlafli's, [31], concerns the variation of the volume (or 
content) of a spherical polytope, 11^, as the defining parameters, such as the po- 
sitions of the vertices, are altered slightly. For arbitrary small displacements, the 
change in volume is 

d\U d \ = ^- [ J2\ U d-2\d9 (60) 

where the sum is over the facets, Tld-2, of the various bounding n^-i's. (These 
facets are Schlafli's secondary perischemons.) The angle 9 is that between the 
(d — 1) -planes intersecting in the facet. (In our previous notation $ij = 9 and 
dM.i fl dM.j is a Iid-2-) Equation (60) can be iterated on d. 

We move quickly to d = 3 and choose a 3-simplex for n 3 , i.e. a general spherical 
tetrahedron, for which the result was also given by Richmond [32], [33]. Denoting 
the edges and corresponding dihedral angles by L and 6>, the volume variation is 

dv = \ E Ldd - ( 61 ) 

edges 

The tetrahedron is uniquely fixed by specifying the six 6>'s or, equivalently, the six 
L's. 

Schlafli points out the following 'reciprocity'. Consider another tetrahedron 
whose angles are n — L and whose edges are n — 9. The change in its volume is 

dV' = -\ J2^~°) dL - ( 62 ) 

edges 

and so d(V + V) is a perfect differential which integrates to 

V + V' = ^-\^-9)L. (63) 

The constant of integration can be determined by setting all the 0's, and L's, equal 
to tv/2 giving the self- reciprocal orthant of volume 7r 2 /8. Schlafli, [34] p 281, fixes 
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the constant by choosing the sides of V vanishingly small, in which case V becomes 
half a full sphere. 

We wish to apply (63) to a doubly rectangular tetrahedron (an orthoscheme), 
such as a fundamental domain, for which the angles are a, (3, 7 and three 7r/2's. 
Schlafli denotes the volume by n 2 f(a, (3, 7)/8, so normalising / to unity for an 
orthant. 

The reciprocal of a doubly rectangular tetrahedron is a tetrahedron (not doubly 
rectangular in general) in which three consecutive 'nonplanar' edges equal tt/2, the 
rest being tt — a, tt — (3 and tt — 7. We write its volume as tt 2 F/8. 

Now we choose a fundamental domain, a = ir/ki, (3 = n/k2, 7 = n/ks, and 
get from (63) 




where on the right we have reverted to our previous notation. Then, from (57), 
using the fact that the volume of 2g fundamental domains is 27r 2 , so that gf = 8, 
we find 

b 2 = g-l-j. (65) 

We recall that, in the two-dimensional (i.e. M 3 ) case, b 2 = g — 1 on the simple 
geometric grounds that rotations fix an axis but not a plane. In M 4 , the elements 
that fix a plane but not a hyperplane will certainly belong to the rotation subgroup 
but now we have the extra possibility of a rotation (a product of four reflections) 
fixing a point but not a plane or hyperplane. These must also be removed, as well 
as the identity. Hence we have 

h = 9 ~ 1 - b 4 , 

which we already know, (45), and whence the curious relation 

F = b 4 f. (66) 

The volume of the reciprocal tetrahedron is 64 times that of the fundamental or- 
thoscheme. Numerically, 64 is the product of the exponents, rrij. For {3 3 }, f = 2/15 
and 64 = 24. At the moment, we have no direct geometrical interpretation of (66). 
The reciprocal orthoscheme does not seem to be a particularly useful construct. 

Schlafli [34] interestingly uses (63), applied to a polytope, to deduce several 
numerical results for /. 

We remark that the Schlafli function, /, has been used occasionally in physics, 
[35], [36], [37]. 
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9. A further identity 

A further identity follows on setting m = 3 in (32), 



t3/2 " r(<i/2 - i) " 1 



_3/2 3/2 

= ^|S d " 3 | ^((d - l)/2 | d) = W |S d " 3 | (26 1 (6 1 - b 2 + d - 3)), (67) 

(see the appendix). This will come into play for d > 3. 

We can only state the result for the (Neumann) coefficient C 3 / 2 in the incom- 
plete form, 

C 3/2 = V f (6tr(^ 2 ) - 3k 2 -4R+ 12(8f - + corner terms (68) 

192 ^JdMi 

where R is the intrinsic curvature of the boundary parts. 
In d-dimensions, the curvature part of C 3 / 2 reduces to 



2 ^{d-l){d-2)\dM\. (69) 

This combines with the last term in (67), if we use (38) for b±, to give for the 
(unknown) corner terms, the rather simple value 

3 /2 

corner terms = - — |S d " 3 | b^h - b 2 - 2) < 0, (70) 
12(7 

in this particular geometry. We check that the right-hand side vanishes for d < 3. 

To the author's knowledge, no general form for the corner term is known. An 
expression in the case that M. is a polygonal cylinder is given by Pathria [38] and 
Baltes, [39], but this is not enough. There are some general results of Cheeger 
[40] . For the time being, we can just take the attitude that we are calculating these 
corner terms for a special geometry. 

The derivation of (40) can be traced back to work of Sommerfeld on propagation 
in a wedge of arbitrary angle. To get the missing terms in (68) it would be necessary 
to discuss a general trihedral corner. 

It would be interesting to use the result (70) to get a handle on the form of 
the corner terms. If d = 3, the right-hand side of (70) must have a combinatorial- 
geometric significance in terms of the vertices of the simplicial decomposition of S 3 . 
In terms of the order g, we find the value of b\{b 2 — b\ + 2) to be 9g/2 for {3 3 }, 6g 



19 



for {3 2 4}, 7g for {343} and 9g for {3 2 5}. For the record, the corresponding values 
of b 2 - h + 2 are 27, 72, 168 and 1080. 

Let us recapitulate the identities. So far we have (36), (37), (47) and (67) (with 
(68)). Assume that the volume of At, and in fact all the geometry of A4, is known 
in terms the dihedral angles. Then, with Schlafli, the order, 2g, of the symmetry 
group can be found from (36). Next, bi follows from (37), and b 2 from (47) but 63 
is curiously absent from (67). 

11. Higher relations 

We return to two dimensions and continue with the discussion in section 4. 
More generally than looking at C(0), we have from the equality of gC,(—k), 



12 1 
Y, mn mCm{-k)-2{Y j n m -l)Cm{-^k-l,-) = - - - - B™ k ( - I did 2 ). 

(71) 



The Bernoulli functions are 'non-topological'. 

We repeat the expression for ( m (—k) given in [12] 

1 2k — 1 m ~ 1 

Cm(-k) = - B 2k+2 (l/2) + ^-^ ( 2 P + l)52fc+i((2p + 1) /2m) . (72) 

^ ' p=0 

The following identity, derived in the appendix to [12], is needed in order to 
rewrite (72) in a form more suitable for substituting into (71). 

k , , x m-1 

ELI (m 2l -l)B 2k _ 2l {-) B 2l = km 2k ~ 2 £ (2p+l) B 2k _ x {(2p+l)/2m) . (73) 
1=0 ^ ' p=0 

For (71) we then find 

>(2) /I 



(2 + fc) 

fc+i 

2k + 



B% k (o\did 2 ) 



Vl E ( 2/C 2 l 2 ) K " l)*«B 2fc+2 _ 2l (i) - 5 2fc+2 (^). (74) 



For k = this is our previous result of course. We write out the explicit form when 
k= 1, 

240Sf } (- I did 2 ) = -4 E n m (2(m 4 - 1) + 5(m 2 - 1)) - 21< 0. (75) 
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Generally, looking at the form of (74), we see that by a recursion argument we 
can derive an expression for the sum Y^i n m( m — 1) in terms of the degrees. Then, 
making use of the relations (3), and of (24), this expression will yield an equation 
for the sums of all odd powers of the three orders, ^ l m 2l ~ 1 , in terms of the two 
degrees. 

There is no doubt that these relations are numerically satisfied. After all they 
are identities. They have nevertheless been checked for the regular solids. (The 
dihedral case [m] may also be used as an example. For this, (p, q, r) = (m, m, 1) 
and di = m, d 2 = 1. Heuristically we can set n p = n m = 1/2 rather than make this 
a special case.) 

-^3 (I I ^1^2) can be calculated or looked up. It comes out directly as 
240^ 2) (i I did 2 ) = -8(4 + di) + 40djdl - UOd^dx + d 2 ) + 

60(4 + 4 + 3did 2 ) - 60(di + d 2 ) + 15 (76) 
but it should be rewritten in terms of g and 61. We have the easily verified relation 

d\ + d\ = b\ + 46? - Agb\ + 6bj + 2g 2 - 8gh - 4g + 4&1 + 1 

and we find for the left-hand side of (75), 

-8bj - 326f + 126? + 286 x + 32#6? - 56#6i + 24# 2 - 28# + 7. (77) 

The resulting identity cannot be obtained from the C-F relation. So far as 
I can see it does not have a geometrical interpretation. Since there are only a 
well defined set of orders (p,q,r) it is perhaps not surprising that there are many 
particular relations. 

12. Conclusion 

Since the approach depends on the expansion of the heat-kernel, it is restricted by 
the availability of explicit forms for the coefficients. It is certainly of interest to 
obtain the unknown corner terms. 

Because of our inability to derive (57) directly (it may of course be obvious), 
we would like to claim a modest amount of novelty for our method. 

At one level, it might be said that all we are doing is checking known coefficients 
in a specific situation. There is, however, a little more to it than that. To check the 
numbers, it is not necessary to rewrite the right-hand side of (40) in terms of the b r . 
It could have been left in the degrees. However the dramatic simplification of the 
Bernoulli functions and the various cancellations suggest some interesting general 
relations that have yet to be uncovered. 



21 



If it is wished to check numerically the identity for d = 4, say, then it will be 
necessary to work out the areas of the 2-faces between neighbouring 3-face cells of 
the boundary of the 4-simplex fundamental domain. 

It does not seem possible to extend easily the methods of this paper to the 
stellated polytopes, such as the Kepler-Poinsot polyhedra. The difficulty is the 
evaluation of the heat-kernels on the branched coverings. In such highly symmetrical 
situations it might be possible to obtain closed forms for the heat-kernels but we 
have not been able to make any progress in this direction. Laporte [20] has some 
useful comments on the mode problem. 

In the two-dimensional case, Rubinowicz [41] claims to have a method of in- 
tegrating the wave-equation on an arbitrary Riemann surface, but, in fact, only 
discusses the case of one branch point in detail and says that the general case fol- 
lows upon a process of 'Zusammenstucklung', which, unfortunately, he does not 
exhibit. 

Appendix 

We give here Sommerville's expressions [28] for the edge lengths of the fundamental 
simplex of the (/-dimensional spherical honeycombs. For {3, 3, . . . , w} = {3 d_1 , w}, 

sin 2 La = - ; ) / , , 78 

3 j + 1 1 - {d - i - 1) cos (2tt/w) V ' 

while, for the reciprocal, {w, . . . , 3, 3} = {to, 3 d-1 }, 

o i — i 1 — dcos (2tt/w) , 

sin 2 L i:i = -rr—. — . / ' . 79 

3 d-i + 1 1 - (j - 1) cos (2tt/w) ' 

For the more general case, {v, 3, 3, . . . , 3, 3, w} = {v, 3 d-2 , w}, 

cqo 2 L _ 1 + (1 + j ~ d) cos (2tv/w) 1 + (1 - i) cos {2ti/v) 

C ° S 13 1+ (1 + i - d) cos (2ir/w) ' 1 + (1 - j) cos (2ir/v) ' ^ ' 

Finally for the case {3, 4, 3, 3, . . . , 3} = {3, 4, 3 d-2 } we have 

9 r d — 7 + 15 — % 

cos 2 L 0j = d ; J + 1 (81) 
5 - J 

and for its reciprocal, {3, 3, ... 3, 4, 3} = {3 d-2 , 4, 3}, 

2 t 5 + j- di + l .. . 

cos^ La = • (i < i < d - 1) 

3 5 + % - d j + 1 v J ~ ' 
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We now turn to the connection of the generalised Bernoulli numbers with the 
Todd Polynomials. 

Hirzebruch [42] gives the relation 

T k (c u ...,c k ) = t^B i k d \d 1 ,...,d d ) k<d (83) 

where, the c s are the elementary symmetric functions of the degrees di (i = 1, . . . , d). 

Using the known formulae for the T k we can obtain the Bernoulli quantities if 
some transformations are made. We wish to write the quantities in terms of the 
elementary symmetric functions, b r , of the exponents rrii (i = 1, . . . , d + 1) so we 
firstly need the connection between the c s and the b r . This is easily found as follows. 
First divide the polynomial, 



d d+1 d 





where 

K = (-!) r+1 bi = (-ir i>i)* b * 

r+1 

using (44). Then set t — > t/(l + t) to get 

d d d 



Expanding the (1+t), the c s can be read off and a rearrangement of the summations 
gives 

c s =b-i)^b-i)^:'I ( m ) 



r=0 



As a point of more than technical interest, Lemmas 1 and 2 of Todd [43] provide 
an alternative means of obtaining the elementary symmetric functions of quantities 
translated by a constant. 

More generally than is needed here we have the following result 

c s {{di +x})=i2^ i ( d+i ~ 
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which follows directly upon the replacement t — > t/(l + Xt) in 

1 o 

or from a Taylor expansion and then use of the standard formula 

dc s ({di}) 



ddi 



sc. 



({*}) 



where the hat signifies omission of the term, (cf Lemma 2 in [43]). Note that Todd's 
h s is a sum of all homogeneous products i.e. 



Eo(-i) s ^ s n(i-**)' 



The quantity we actually require is the Bernoulli function B r (x\d) (at x = 
(d — l)/2) which is a polynomial in x whose coefficients are the Bernoulli numbers 
in (83), 

Bf° {x I d) = (-l)V! £(-1)' ^T r _, (d, . . . , c r _ ; ). (85) 

z=o 

From the listed expressions for the Todd polynomials in the c s , using (84) and 
(85), we can generate the Bernoulli polynomials at x = (d — l)/2 in terms of the b r . 
These have been used in the main text. We collect a few results here. 

B { 2 d) ((d - l)/2 I d) = | {bi(h - 1) + b 2 + (3 - d)/2) 

B ( 3 d) ((d -l)/2 | d) = hb! (h -b 2 + (d- 3)/2) (86) 

B{ d) ((d - l)/2 | d) = — X — (86 4 - 246 3 - 863&1 - 246^ - 46 2 (86? - I661 - 5d + 9)+ 

46i(6i - 1)(26? + IO61 + 5rf - 13) - (5d - 3)(d - 5)). 

It is probably more instructive to evaluate the Bernoulli polynomials directly 
rather than go through the Todd polynomials. We can use the general theory 
of multiplicative series as developed by Hirzebruch [42] and this will be detailed 
elsewhere. 
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